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Abstract 

In recent years, there has been considerable interest in understanding the 
motion in Hamiltonian systems when phase space is divided into stochastic 
and integrable regions. This paper studies one aspect of this problem, namely, 
the motion of trajectories in the stochastic sea when there is a small island 
present. The results show that the particle can be stuck close to the island 
for very long times. For the standard mapping, where accelerator modes 
are possible, it appears that the mean squared displacement of particles in 
the stochastic sea may increase faster than linearly with time indicating non- 
diffusive behavior. 



Introduction 

Many important problems in physics are described by Hamiltonians of two 
degrees of freedom. Examples are the motion of a charged particle in electro- 
static waves, the motion of a charged particle in various magnetic confinement 
devices, the acceleration of a particle bouncing between a fixed and an oscil- 
lating wall, the wandering of magnetic field lines, etc. In such systems, there 
is usually a range of parameters (normally when the coupling between the two 
degrees of freedom is large), where the motion in nearly the whole of phase 
space is stochastic. Such behavior is seen for instance in the standard mapping 
(Chirikov, 1979), 

r t -r t -i = -ksin.9 t _i, 6 t -9 t _ 1 = r t . 

When k is large most of phase space is stochastic. However, there may still be 
small islands present; stochastic trajectories can wander close to these islands 
and remain there for a long time leading to unexpectedly long correlations. 
The effect of these correlations can be dramatic. The simplest approximation 
for the diffusion coefficient, 
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(where the average is over some appropriate ensemble), is given by assuming 
that the phase 9 is a random variable in the equation for r. This gives the 
"quasi-linear" result V = V q i = \k 2 . However, a numerical determination 
(Karney et al, 1982) of the diffusion at k = 6.6, where the random phase 
approximation might be expected to be accurate, gave V/Vqi ~ 80. At this 
value of k there is an island ("accelerator mode") present in the stochastic sea. 
This leads to long-time correlations in the acceleration of the particle and an 
enhanced diffusion coefficient. 

In this paper, we examine more closely the effect these islands have on 
a stochastic trajectory. As far as determining the effect on the correlation 
function, this involves determining how "sticky" the island is. Given that the 
stochastic trajectory comes within a certain distance of the boundary of the 
island, how long do we expect it to stay close to the island? This approach 
is inspired by work of Channon and Lebowitz (1981) on the correlations of a 
trajectory in the stochastic band trapped between two KAM surfaces in the 
Henon map. Similar work has been carried out on the whisker map by Chirikov 
and Shepelyansky (1981); this work is being extended by B. V. Chirikov and 
F. Vivaldi. The work reported herein is described in more detail by Karney 
(1983). 

Derivation of Mapping 

Far into the stochastic regime for a general mapping, the islands which appear 
via tangent bifurcations are very small and exist only for a small interval in 
parameter space. This allows us to approximate them by a Taylor expansion in 
both phase and parameter space about the tangent bifurcation point retaining 
only the leading terms. This was carried out by Karney et al. (1982) where 
the resulting mapping was reduced to a canonical form 

Q ■ Ut — Vt-i = 2(xj_ 1 — K), x t -x t -i = y t . 

Here K is proportional to k — /c tang (fetang is the parameter value where the 
tangent bifurcation takes place) and x and y are related to the original phase 
space coordinates by a smooth transformation. The quadratic mapping Q rep- 
resents an approximation of the general mapping close to the point of tangent 
bifurcation. For < K < 1, this mapping has stable (elliptic) and unstable 
(hyperbolic) fixed points at (x,y) = (^f\/ r K,0), respectively. The elliptic fixed 
point is usually surrounded by integrable trajectories (KAM curves) which de- 
fine a stable region (the island) in which the motion is bounded. An example 
of island structure is shown in Fig. 1 for K = 0.1 (the value of K at which 
extensive numerical calculations have been carried out). 

Referring to the islands shown in Fig. 1, consider a particle which at t = 
is close to, but outside, the islands. Initially, the particle will stay close to 
the islands; however as we let t — > ±oo, we find (x,y) — > (oo, ±oo). It is just 
such trajectories we are interested in, because they correspond to particles in 
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Fig. 1. Some islands of the quadratic map Q for K = 0.1. 

the stochastic region of the general mapping approaching the islands, staying 
there for some time (and contributing to long-time correlations), and then 
escaping back to the main part of the stochastic region. 

What we need is some way of bringing these particles back to the vicinity 
of the island. We do this by defining an L x L square around the island. 
This square spans the region x m - in < x < x max and —\L<y<\L where 
L = £ m ax — x m i n . Whenever an orbit leaves this square at (x t ,y t ), we pick a 
new initial condition 

(x , yo) = {x t - mL, y t - nL) 

with m and n being integers chosen so that (x , y ) lies inside the square. We 
also record the length t of the previous orbit segment. This procedure defines 
the periodic quadratic map Q*, which can be shown to sample the orbits close 
to the island in the same way that the general mapping does. Examples of the 
orbits of Q* are shown in Fig. 2 for the same parameters as for Fig. 1. 

One useful way of looking at Q* is as a magnification of a small region 
near a tangent bifurcation in the general mapping. The difference is that once 
the trajectory leaves the vicinity of the islands, it is immediately re-injected 
on the other side of the islands. In the general map, the trajectory will spend 
some long time, which depends on the ratio of the size of the islands to the 
total accessible portion of phase space, in the stochastic sea before coming 
back to the vicinity of the islands. 

Assuming that the long-time behavior of stochastic orbits is dominated 
by the region close to the islands, there are two advantages to reducing the 
problem to a study of Q*. Firstly, since Q* describes the behavior of most 
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Fig. 2. (a) Stochastic trajectories for periodic quadratic map Q* for K = 0.1. 
(b) An enlargement of a portion of (a). Here x m i n + \^K = —0.4, x max + \f~K = 
0.6. 

islands far into the stochastic regime, the properties of many mappings may 
be treated by looking at a special mapping Q* which depends only on a single 
parameter K. The second advantage is that the properties of orbits close to 
the islands may be studied much more efficiently because there is no need 
to follow orbits while they spend a long and uninteresting time far from the 
islands. 
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Trapping Statistics 

The prescription for numerically determining the stickiness of the island sys- 
tem in Q is to compute a long orbit in the stochastic region of Q*. The orbit 
is divided into segments at those points where it leaves the L x L square. The 
main results of the calculation are then the trapping statistics f t which are 
proportional to the number of orbit segments which have a length of t. Sup- 
pose that the total length of the orbit is T and N t is the number of segments of 
length t. If T is so large that we can ignore partial segments at the ends of the 
orbit, then we have ^tN t = T; the total number of segments is N = Y2^t- 
The trapping statistics are defined by ft = N t /T and are therefore normalized 
so that Yl^ft = 1- The mean length of the orbits is given by a = l/X)/t 
(= T/N). The probability that a particular segment has length t is p t = aft 
(= Nt/N). If an arbitrary point is chosen in the orbit, then tft is the proba- 
bility that this point belongs to a segment of length t and ft is the probability 
that it belongs to the beginning, say, of a segment of length@t. 
The survival probability 

oo 

Pt= £ Pr 

T = t+1 

is the probability that an orbit beginning in a segment at t = is still trapped 
in the same segment at time t. Note that Pq = 1 as required. This is the quan- 
tity studied by Channon and Lebowitz (1980) and Chirikov and Shepelyansky 
(1981). The correlation function 

oo oo 

C T = J2(t-T)ft = J2Pt/a 

t=T t = T 

is the probability that a particle is trapped in the same segment at two times 
r apart. Again, we have Co = 1. 

There is another way of interpreting C T : Consider a drunkard who executes 
a one-dimensional random walk with velocity v = dr/dt = ±1. The direction 
of each step is chosen randomly, while the durations of the steps are chosen to 
be the lengths of consecutive trapped segments of Q* . Then for integer r, C T 
is just the usual correlation function for such a process, i.e., [vtVt+ T )t- The 
behavior of this random-walk process is similar to the behavior of an orbit 
in the general mapping when two accelerator modes with opposite values of 
the acceleration are present. (This is the case with the first-order accelerator 
modes for the standard mapping.) 

A diffusion coefficient may be defined by 

oo 

T = l 

This gives the diffusion rate for the drunkard in the random-walk problem 
above. It is also related to the diffusion coefficient for the general mapping. 
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Results for K = 0.1 

We have measured f t for K between and 1.3 at intervals of 0.05, and at most 
of the values of K a slow algebraic decay of f t is seen. A representative case 
is K = 0.1, whose trapping statistics are given in Fig. 3(a), which illustrates 
the slow decay for very long times t ~ 10 7 . Also given in Fig. 3 are Pt, C T , 
and a = — d log C T /d log r (thus locally C T ~ T~ a ). This last plot shows the 
power at which C T decays varying between about j and |. 

A glance at Fig. 2 shows the origin of this behavior. The central island is 
surrounded by a chain of sixth-order islands. Around each of these islands are 
several other sets of islands. This picture repeats itself at deeper and deeper 
levels. A particle which manages to penetrate into this maze can get stuck in 
it for a long time. 

For r < 10 4 , Fig. 3(d) gives a « |. Correspondingly we have P t ~ t~ p 
where p = 1 + a w 5/4. This is close to the asymptotic (i — ► oo) result 
found by Chirikov and Shepelyansky (1981) for the whisker map, in which 
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(p) ps 1.45. However, in our case, a shows some strong variations beyond 
r ps 10 4 where C T "steps down" (e.g., between 10 4 and 3 x 10 5 ). This means 
that the asymptotic form of C T is very difficult to determine numerically. 

The diffusion coefficient D is given by the summation of C and is approx- 
imately 6400. The error in this estimate of D depends on the asymptotic 
form for C T . On the basis of the numerical results, we cannot rule out the 
possibility that as r — > oo, C T decays with a < 1. In that case, D would be 
infinite! 

If D is indeed infinite, we would wish to know how a group of particles 
spreads with time. We again consider the drunkard's walk based on Q* which 
was introduced earlier. The second moment of r is related to the correlation 
function by 

t 

St = <(r t - r ) 2 ) = tC + 2 £(t - r)C T . 

T=l 

This is plotted in Fig. 4(a), using the data of Fig. 3. For t < 10 4 , St grows 
somewhat faster than i 3 / 2 (see Fig. 4(b)) and even until t ps 10 7 , S t is growing 
significantly faster than linearly. Beyond 10 7 , the numerical data shows a 
convergence to a linear rate; but this is merely because no segments longer 
than about 6 x 10 7 were observed. For t — > oo, St grows as t 2 ~ a , assuming 
that the exponent a at which C T decays asymptotically is less than 1. If the 
diffusion coefficient is estimated from D t = ^St/t, then D t grows with t as 
shown in Fig. 4(c). 

Discussion 

We can apply these results to the determination of the correlation function of 
a general mapping. Suppose the correlation function is defined by 

C(r) = (/i(x(t))/i(x(t + r))) t) 

where h is some smooth function of the position in phase space x. Then the 
contribution of an island located at xq to C(t) is (Karney, 1983) 

Cis(r) = /i 2 (x )(B/A)C T , 

where A is the total area of the stochastic component of the general mapping 
and B is the portion of the that area which is near the island. (More precisely, 
when we regard the L x L square of Q* as being a magnification of a small 
area of the general mapping, then B measures the area of the stochastic com- 
ponent within this small area.) Similar relations connect D and St and the 
corresponding quantities for the general mapping. 

In the case of accelerator mode in the standard mapping, K is related to 
the parameter k by k 2 = (2irn) 2 + 16K where n is an integer. For K = 0.1, 
we take 6400 as a lower bound for D. We find that the contribution to the 
diffusion coefficient is increased over its quasi-linear value by a factor of at 
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least 360/n 2 . Thus for n = 1 or k ~ 6.41, the islands completely dominate 
the diffusion. The first-order accelerator modes continue to have such a large 
effect at least until k ~ 100. If D is in fact infinite, even arbitrarily small 
accelerator modes will eventually dominate the motion and Fig. 4 can be used 
to estimate the time at which the accelerator modes become important. 

In summary, small islands within the stochastic sea lead to correlations in 
the stochastic orbits for extremely long times. When the islands are accelerator 
modes, this may cause the particles to behave non-diffusively, i.e., the mean 
squared displacement of the particles may increase faster than linearly with 
time. 

In order to provide a definitive answer to this question, the asymptotic 
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behavior of C T must be determined. Because the asymptotic regime starts at 
such a large r (greater than 10 7 ), it appears its properties cannot be studied 
by the numerical method used in this paper. What is needed is a better 
analytical understanding of the behavior of trajectories close to the border 
between integrability and stochasticity. Chirikov (1982) has made some useful 
steps in this direction. 
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